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PERIODIC JITTER CHARACTERIZATION USING PSEUDO- 
RANDOM SAMPLING 

Background of the Invention 

High speed signals in communication systems are often corrupted by signal timing 
instabilities, typically referred to as signal jitter. Signal jitter can include various 
components, such as random jitter, data dependent jitter, and periodic jitter. The periodic 
jitter can be characterized by determining the frequency of the periodic jitter, which is 
useful for identifying the source of the periodic jitter. 

Sampling systems, such as equivalent-time sampling oscilloscopes, are well-suited 
for characterizing high speed signals that have repetitive bit patterns. Highly stable time 
bases within the sampling systems also make the sampling systems suitable for 
characterizing periodic jitter. However, the sampling rates of the sampling systems are 
typically lower than the bit rates of the high speed signals and the frequencies of the 
periodic jitter that may be present on the signals. This under-sampling results in aliasing, 
which impairs the ability of the sampling system to characterize periodic jitter. 
Particularly, this aliasing makes it difficult for the sampling system to distinguish between 
actual signal components of the periodic jitter and aliased signal components that result 
from the under-sampling. Accordingly, there is a need for a scheme for characterizing 
periodic jitter that distinguishes actual signal components of periodic jitter from the aliased 
signal components of the periodic jitter. 
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Summary of the Invention 

A system and method according to the embodiments of the present invention 
characterize jitter of an applied signal. The characterization includes acquiring a set of 
pseudo-randomly timed samples at a designated position on the signal, assigning a jitter 
5 value to each of the pseudo-randomly timed samples in the acquired set, and selecting a 
frequency from an array of frequencies based on a correlation of the assigned jitter values 
with the frequencies in the array. The periodic jitter associated with the signal is 
designated to have the frequency within the array of frequencies that has the highest 
correlation to the assigned jitter values. 

10 

Brief Description of the Drawings 

Figure 1 shows a sampling system for characterizing periodic jitter of repetitive 
signals according to embodiments of the present invention. 

Figure 2 shows a method for characterizing periodic jitter according to 
1 5 embodiments of the present invention. 

Figure 3 shows a detailed view of an alternative embodiment of the method of 
Figure 2. 

Figures 4A-4F show exemplary waveforms at various stages in the 
characterization of periodic jitter according to embodiments of the present invention. 

20 

Detailed Description of the Embodiments 

Figure 1 shows a sampling system 10 suitable for characterizing periodic jitter of a 
signal SI according to embodiments of the present invention. The signal SI is typically a 
high speed data signal, or other signal, having an associated bit rate Rl and pattern length 
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PL While the signal SI is either repetitive, or semi-repetitive with one or more sequences 

that repeat with each pattern length PI, the signal SI is hereinafter referred to as "repetitive 
signal SI". A source 12, typically a communication system, network node, instrument, 
device or other element, provides the repetitive signal SI to a sampler 14 within the 
5 sampling system 10. A pattern trigger Tl, typically provided by the source 12 or derived 
from the repetitive signal SI, is applied to a trigger input TRIG of the sampler 14, 
providing a timing reference for the acquisition of samples by the sampler 14. 

The sampler 14 is an equivalent-time sampling oscilloscope, signal digitizer, 
analog-to-digital converter or other suitable sample acquisition system. The sampler 14 

1 o has a sampling rate R s that is typically lower than the bit rate Rl of the repetitive signal S 1 
and the frequencies of the periodic jitter that may be present on the repetitive signal SL 
For example, when the sampler 14 is an AGILENT TECHNOLOGIES, INC. model 
86100 DCA, a typical sampling rate R s is 40 kHz, whereas the bit rate Rl of a typical 
repetitive signal SI may be as high as 10 Gbits/second and the periodic jitter may have 

15 frequencies as high as 100 MHz. The sampling system 10 includes a controller 16 that 
initiates sample acquisitions according to the pattern trigger Tl and processes acquired 
samples to characterize the periodic jitter present on the repetitive signal SL 

Figure 2 shows a method 20 for characterizing periodic jitter, suitably 
implemented with the sampling system 10. The method 20 includes establishing an array F 

20 of frequencies f (step 22), acquiring a set S r of pseudo-randomly timed samples at a 
designated position 17 in the repetitive signal SI (step 24), assigning a jitter value to each 
of the samples in the acquired set S r (step 26) and selecting a frequency f j from the array F 
of frequencies f based on correlation of the assigned jitter values with frequencies f in the 
array F (step 28). The periodic jitter associated with the repetitive signal SI is designated 
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to have the frequency fj, which is the frequency within the array F of frequencies f that has 

the highest correlation. 

According to one embodiment of the present invention, the array F of frequencies f 
in step 22 includes frequencies identified as suspect, for example, frequencies of 
5 interference signals or interference sources, oscillation frequencies of components within 
the source 12 providing the repetitive signal SI, or other identified frequencies of periodic 
jitter imposed on the repetitive signal SI. In an alternative embodiment of the present 
invention, the array F of frequencies f includes the frequency fpEAK of a signal peak S PE ak 
in a spectrum SPER (shown in Figure 4A) resulting from uniform periodic sampling at a 

10 designated position 19 of the repetitive signal SI at a sampling rate Rs that is lower than 
the bit rate Rl of the repetitive signal SI, and frequencies f that are offset from the 
frequency fpEAK by integer multiples n of the sampling rate R s . Thus, in this embodiment, 
the array F includes the frequencies {fpEAK ± n R s }, where the integer multiple n has 
integer values such that 0 < n < R1/(2RS). 

1 5 Figure 3 shows exemplary steps 22a-22e for establishing the array F of frequencies 

f in step 22 in the embodiment where the array F includes the set of frequencies {fpEAK ± n 
Rx} based on the signal peak S PE ak in the spectrum SPER. Steps 22a and 22b comprise 
characterizing the relationship between amplitude and time on a designated amplitude 
transition 13 of the repetitive signal SI. Typically, the amplitude transition 13 is a rising or 

20 falling edge transition between logic states encoded in the bit stream of the repetitive signal 
SI. However, the amplitude transition 13 is alternatively any suitable feature of the 
repetitive signal SI exhibiting an amplitude change versus time that repeats with the 
pattern length PL 
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In step 22a, the designated amplitude transition 13 is characterized. This 

characterization includes, for example, acquiring samples {Si...S N } at various times 
{ti...t N } along a repeating edge transition in the repetitive signal SI according to the 
pattern trigger Tl, using equivalent-time sampling techniques. Figure 4B shows an 
5 example of samples {Si . . .S N } acquired along a rising edge transition. 

In step 22b, a mapping between amplitudes {Ai...A N } of the acquired samples 
{Si...S N } and the corresponding times { ti — t N } of the acquired samples {Si...Sn} is 
established. Typically, this mapping involves a linear function relating the amplitudes 
{Ai...A N } and times {ti...t N }, a polynomial relating the amplitudes {Ai...A N } and times 
10 {ti...t N }, or a look-up table relating the amplitudes {Ai...A N } and times {ti...t N }. 
However, any number of suitable techniques are alternatively used to establish this 
mapping between the amplitudes {Ai...A N } and times {ti...t N } of the acquired samples 
{Si...Sn}. 

In step 22c, a set Su of samples {Syi... Sux} is acquired at a designated position 
15 19 on the amplitude transition 13 characterized in step 22a. Consistent with the example 
above, the designated position 19 is the midpoint of the rising edge transition shown in 
Figure 4B. The designated position 19 is at a nominal time position tu NO M that has a 
corresponding nominal amplitude Aunom- Timing instabilities of the repetitive signal SI, 
including periodic jitter, result in the amplitudes { Aui . . . Aux} of acquired samples {Sui . . . 
20 Sux} deviating from the nominal amplitude Aunom, as shown in Figures 4C-4D. Figure 
4C is an equivalent-time representation of the repetitive signal SI showing repeating 
occurrences of the amplitude transition 13 on which the designated position 19 is located. 

While the timing instabilities vary the timing of the repetitive signal SI, the highly 
stable time base of the sampler 14 within the sampling system 10 enables the acquired 
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samples {Sui... Sux} within the set Su to be periodically acquired at times {tui...tu x } 
spaced by precisely defined, uniform time intervals tu- The time intervals tu between the 
times {tui...tu x } are each equal to a fixed integer multiple M of the pattern length PI 
associated with the repetitive signal SI. Thus, samples {Sui... Sux} within the set Su are 
5 acquired every Mth repetition of the pattern length PI as shown in the representation of the 
acquired samples {Sui... Sux} of Figure 4D. 

In step 22d, a corresponding jitter value is assigned to each of the samples {Sui... 
Sux} in the set Su of periodically acquired samples in step 22c, resulting in a set Ju of jitter 
values. This assignment includes determining for the samples {Sui... Sux} in the set Su, 

10 the deviations {Dui...Dux} of the amplitudes {Aui-..Aux} of each of the samples {Sui... 
Sux} from the nominal amplitude Aunom. The mapping of amplitudes {Ai...A N } and 
times { ti — t N } established in step 22b is then used to convert the determined amplitude 
deviations {Dui - Dux} to corresponding timing deviations. The resulting timing 
deviations comprise the set Ju of jitter values. When a linear mapping is established in 

15 step 22b, the set Ju of jitter values is obtained by dividing the amplitude deviations 
{Dui...Dux} by the slope of the linear function relating amplitudes {Ai...A N } and times 
{tj...t N }. However, when a polynomial mapping is established in step 22b, the set Ju of 
jitter values is obtained by evaluating the polynomial for each of the amplitude deviations 
{Dui . . .Dux} from the nominal amplitude Aunom- When the mapping in step 22b is a look- 

20 up table, the set Ju of jitter values is assigned according to the look-up table, typically 
using interpolation to accommodate deviations from the nominal amplitude Aunom that 
fall between values in the look-up table. 

In step 22e, the set Ju of jitter values resulting from step 22d is transformed to the 
corresponding spectrum SPER using a Discrete Fourier Transform or other suitable 
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transform between the time domain and the frequency domain. An exemplary spectrum 

SPER based on the set Sy of samples {Sui... Sux} is shown in Figure 4A. From this 
spectrum SPER, the array F of frequencies f is established for the embodiment of the 
present invention in which the array F includes the frequency fpEAK of an identified signal 
5 peak Speak within the spectrum SPER and frequencies offset from the frequency fpEAK of 
the signal peak S PE ak by the integer multiples n of the sampling rate R s at which the 
samples {Sui-.. Sux} in step 22c are periodically acquired. While an array F is shown 
associated with one of the signal peaks S PE ak in the spectrum SPER, a separate array F can 
be established based on other signal peaks identified in the spectrum. 

10 Once the array F is established in step 22 via steps 22a-22e or other means 

according to the above-recited embodiments of the present invention, the remaining steps 
24-28 of the method 20 are executed. In step 24, the set S r of pseudo-randomly timed 
samples is acquired at the designated position 17 on a designated amplitude transition of 
the repetitive signal SI as shown in Figure 4E. In the example where the array F of 

15 frequencies f is established according to steps 22a-22e, the designated position 17 can be 
timed by the sampler 14 to lie on the amplitude transition 13 characterized in steps 22a- 
22b so that the relationship between amplitudes {Ai...A N } and times {ti...tN} of the 
acquired samples {Si . . .S N } established in steps 22a and 22b can be used in the assignment 
of jitter values in step 26. When the designated position 17 lies on this amplitude 

20 transition 13, periodic jitter associated with the repetitive signal SI can also distinguished 
from data dependent jitter associated with the repetitive signal SI. 

Figure 4E is an equivalent-time representation of the repetitive signal SI showing 
repeating occurrences of the amplitude transition on which the designated position 17 is 
located. In the exemplary waveform represented in Figure 4E, the designated position 17 
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is shown on the midpoint of a rising edge transition. The designated position 17 is at a 

nominal time position tr NO M on this amplitude transition that has a corresponding nominal 
amplitude Arnom. However, timing instabilities of the repetitive signal SI, including 
periodic jitter, result in the amplitudes {Ari...Ark} of the acquired samples {Sri... Srk} 
5 within the set S r deviating from the nominal amplitude Arnom, as shown in Figures 4E-4F. 
While the timing instabilities vary the timing of the repetitive signal SI, the highly stable 
time base of the sampler 14 within the sampling system 10 enables the acquired samples 
{ Sri . . . Srk} within the set S r to be acquired at precisely defined times t k . 

The samples {Sri... Srk} in the set S r are timed according to the relationship of 
10 equation 1, where t k represents the timing of the kth acquired sample, L is an integer, r k is a 
random integer associated with the kth acquired sample, and where PI is the pattern length 
PI associated with the repetitive signal SI and Rl is the bit rate Rl associated with the 
repetitive signal SI, as previously designated. 

'Pp 



tk + 1 = tk + 



(L + rk) (1) 



1 5 Equation 1 shows that successive samples t k , t k+ i within the set S r are acquired at 

time intervals tr that are pseudo-random integer multiples L+r k of the ratio of the pattern 
length PI to the bit rate Rl. 

In step 26 of the method 20, a jitter value Jitterik) is assigned to each of the 
samples {Sri. . . Srk} in the set S r , where k is an integer such that 0<k<K, to provide the 

20 corresponding set J r of jitter values. This assignment includes determining for the samples 
{Sri... Srk} in the set S r , deviations {Dri...Drk} of the amplitudes {Ari...A rk } of the 
samples from the nominal amplitude Arnom (shown in Figure 4F) and includes converting 
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the determined amplitude deviations {Dri . . .Drk} to corresponding timing deviations. The 

resulting timing deviations comprise the set J r of jitter values Jitter(k). 

This assignment is based on the relationship between amplitude and time on the 
amplitude transition of the repetitive signal SI upon which the designated position 17 is 
5 positioned, and is determined analogously to the assignment of jitter values in step 22d to 
the periodically timed samples acquired in step 22c, with the exception that here, the 
samples {Sri... Srk} are acquired at the non-uniform time intervals tr determined by 
equation 1. When the designated position 17 coincides with the amplitude transition 13 on 
which the mapping of amplitudes { Ai . . . A N } and times {ti . . .t N } of step 22b is established, 

10 this mapping can be used in assigning the jitter values to the set S r of samples {Sri... 
Srk}. For example, when a linear mapping is established in step 22b, the set J r of jitter 
values is obtained by dividing the amplitude deviations {Dri... Drk} by the slope of the 
linear function relating amplitudes {Ai...A N } and times {ti...t N } on the characterized 
amplitude transition 13. When a polynomial mapping is established in step 22b, the set J r 

15 of jitter values is obtained by evaluating the polynomial for each of the deviations 
{Dri... Drk} of the amplitudes {Ari...A rk } from the nominal amplitude Arnom- When 
the mapping in step 22b is a look-up table, the set J r of jitter values is assigned according to 
the look-up table, typically using interpolation to accommodate amplitude deviations 
{Dri... Drk} from the nominal amplitude Arnom that fall between values in the look-up 

20 table as needed. 

However, when the designated position 17 does not coincide with the amplitude 
transition 13 on which the mapping between amplitudes {Ai...A N } and times {ti...t N } is 
established in steps 22a-22b, a separate mapping between amplitude and time can be 
developed for the amplitude transition that contains the designated position 17 in a manner 
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However, when the designated position 17 does not coincide with the amplitude 

transition 13 on which the mapping between amplitudes {Ai...A N } and times {tj...t N } is 
established in steps 22a-22b, a separate mapping between amplitude and time can be 
developed for the amplitude transition that contains the designated position 17 in a manner 
similar to steps 22a and 22b. This mapping can then be used to convert the amplitude 
deviations { A R! . . . Ark} to the timing deviations that comprise the set J r of jitter values. 

In step 28 of the method 20, a frequency is selected from the array F of 
frequencies f based on a correlation of the assigned jitter values Jitterik) in the set J r with 
frequencies f in the array F. An exemplary correlation is shown in equation 2, where K is 
the number of acquired samples {Sri... Srk} in the set S r , and where f represents the 
frequencies in the array F. The periodic jitter associated with the repetitive signal SI is 
designated to have the frequency fj, which is the frequency within the array F of 
frequencies f that has the highest correlation R(f) to the assigned jitter values Jitterik) in 
the set J r .. 



R(f) = 



-^^Jitter(k).e j2 " tkf 



(2) 



When a linear function relates amplitudes and times on the amplitude transition 
that contains the designated position 17, the amplitudes {Ari...Ark} of the samples 
{Sri... Srk} and the jitter values Jitterik) in the set J r of jitter values are related by a 
constant. Based on this relationship, in one embodiment of the present invention, the 
assignment of jitter values to the samples {Sri... Srk} in step 26 comprises determining 
the amplitudes {Ari...Ark} of the acquired samples {Sri... Srk}, and the selection of the 
frequency fj from the array F of frequencies f based on a correlation of the assigned jitter 
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frequencies f that has the highest correlation to the amplitudes of the acquired 

samples{ Sri... S RK }. 

While the embodiments of the present invention have been illustrated in detail, it 
should be apparent that modifications and adaptations to these embodiments may occur to 
5 one skilled in the art without departing from the scope of the present invention as set forth 
in the following claims. 
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